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We study dual superconductivity of the ground state of SU(2) gauge theory, in connection with confinement. We do that 
measuring on the lattice a disorder parameter describing condensation of monopoles. Confinement appears as a transition to 
dual superconductor, independent of the abelian projection defining monopoles. Some speculations are made on the existence 
of a more appropriate disorder parameter. A similar study for SU(3) is presented in a companion paper. 
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I. INTRODUCTION 

Order-disorder duality pl,B] plays an increasingly im- 
portant role in our understanding of the dynamics of 
gauge theories, specifically of QCD fipf and of its su- 
persymmetric generalizations J5J. 

Duality is typical of systems which can have configu- 
rations with non trivial spatial topology, carrying a con- 
served topological charge. The prototype example is the 
2d Ising model. If viewed as a discretized version of a 1 + 1 
dimensional field theory, it presents one-dimensional con- 
figurations, kinks, whose topology is determined by the 
boundary conditions (±1) at x\ = ±oo. 

In the usual description in terms of the local variable 
ct(.t) = ±1, at low temperature (weak coupling) the sys- 
tem is in an ordered phase with nonzero magnetization 
(a) ^ 0. At the critical point (a) — » and the system 
becomes disordered, (a) is called an order parameter. 
However one can describe the system in terms of a dual 
variable a*, on a dual lattice. A dual Id configuration 
with one a* up is a kink, which is a highly non local 
object in terms of a. In ref. H it was shown that the 
partition function in terms of a* has the same form as 
in terms of a, i.e. that the system with dual descrip- 
tion looks again as an Ising model, except that the new 
Boltzmann factor (3* is related to the old one /3 by the 
relation 



sinh (2/3*) 



sinh (2/3) 



(1) 



The disordered phase is an ordered phase for the dual 
and viceversa. In the disordered phase (a*) ^ 0: kinks 
condense in the ground state, (a*) is called a disorder 
parameter, a* is a dual variable to a. In this specific 
case the system is self dual, and duality transformation 
maps the strong coupling regime in the weak coupling 
regime and viceversa. 



Other systems showing duality properties are the 3d XY 
model, whose dual is a Coulomb gas in 3d, and the com- 
pact U(l) gauge theory. 

In the 3d XY model topological excitations are the 
vortices of the 2d XY model. These vortices condense in 
the disordered phase g . 

For U(l) theory topological excitations are monopoles. 
There the duality transformation can be performed for 
special choices of the action (e.g. the Villain action [R, 
dual to a Z gauge theory, and the Wilson action g). For 
other choices it is not known how to explicitly perform 
the transformation to the dual. 

An alternative approach consists in identifying the 
symmetry which is spontaneously broken in the disor- 
dered phase, i.e. the topological configurations which are 
supposed to condense, and in writing a disorder param- 
eter in terms of the original local fields || . The disorder 
parameter is then the vacuum expectation value (vev) of 
a non local operator. 

This approach has been translated on the lattice 
jniinj, tested by numerical simulations in the compact 
U(l) gauge theory |[12| , in the 3d XY model || and in 
the 0(3) sigma model |13|] , and first used to investigate 
colour confinement in QCD in ref. [ [l4| . 

In the early literature on the subject condensation was 
demonstrated as the sudden increase of the density of 
topological excitations. This is incorrect, since disorder 
can only be described by the vev of an operator which 
violates the dual symmetry and the number of excitations 
does not. 

Looking at symmetry is specially important in QCD. 
For QCD there exists some general idea about the dual 
P)[l5|. The dual description should also be a gauge the- 
ory, possibly with interchange of the role of electric and 
magnetic quantities. 

This idea could fit the mechanism for confinement of 
colour proposed in ref.'s [|l6|,[l7f as dual superconductiv- 
ity of the ground state, if confinement were due to disor- 
der and monopoles were the topological excitations which 
condense. However a dual superconductor is a typically 



abelian system, while the disorder parameter is expected 
to break a non abelian symmetry. An abelian conserved 
monopole charge can be associated to each operator in 
the adjoint representation by a procedure which is known 
as abelian projection Jig] . We will recall that proce- 
dure in Sect. |l|. There exists a functional infinity of 
choices for the operator, and correspondingly an infin- 
ity of monopole species. A possibility is that the true 
disorder symmetry implies the condensation of all these 
species of monopoles [[L8[ . Some people believe instead 
that some abelian projection (specifically the Maximal 
Abelian one) identifies monopoles that are more relevant 
than others for confinement. Both attitudes reflect our 
ignorance of the dual description of the system. 

In this paper we shall systematically explore conden- 
sation of monopoles defined by different abelian projec- 
tions, in connection with confinement of colour. 

We will do that for SU(2) gauge theory. The treatment 
of SU(3) will be given in a companion paper. Some of the 
results have been obtained during the last years and have 
been reported to conferences and workshops ]iq , |l9| , p0{ . 
This paper contains conclusive results, and is an organic 
report of the methods and of the results obtained after 
ref. Q. 

Our strategy consists in constructing an operator with 
non z ero magnetic charge, for each abelian projection 
(Sect. III). Its vev is a candidate disorder parameter for 
dual superconductivity of the ground state. We shall de- 
termine numerically that vev at finite temperature below 
and above the deconfining phase transition. If condensa- 
tion of these monopoles is related to confinement, we ex- 
pect the disorder parameter to be zero in the deconfincd 
phase, and different from zero in the confined phase. 

This is strictly speaking true only in the thermody- 
namic (infinite volume) limit. A finite size scaling analy- 
sis allows to go to that limit, and, as a by-product, gives 
a determination of the transition temperature and of the 
critical indices if the transition is higher order than first. 
This analysis is presented in Sect. IV. 

A special treatment for the Maximal Abelian projec- 
tion is presented in Sect. [V[ 

We find that gauge theory vacuum is indeed a dual su- 
perconductor in the confined phase, and becomes normal 
in the deconfined phase for a number of abelian projec- 
tions, actually for all projections that we have analyzed, 
in agreement with the guess of ref. pll . 

The idea that confinement is produced by dual super- 
conductivity is thus definitely confirmed. The guess that 
all the abelian projections are physically equivalent is also 
supported, and this is an important piece of information 
on the way to understand the true dual symmetry. 

We find evidence that 517(2) deconfining transition is 
second order. In next paper we will show that for SU(3) 
this transition is first order. 

An analysis of full QCD, including quarks, is on the 
way; if the mechanism proved to be the same, the idea 
that quarks are a kind of perturbation, and that the dy- 
namics is determined by gluons would be tested. This 



would also be a test of the ansatz that the theory al- 
ready contains its essential dynamics at N c — oo, and 
that the presence of fermions and the extrapolation to 
N c = 3 can be viewed as perturbations. 
The results are summarized in Sect. N% 



II. THE ABELIAN PROJECTION 

What follows will refer to the case of gauge group 
SU(2). Adaptation to £{7(3) will be described in the 
companion paper. 

Let $(a;) be the direction in colour space of any local 
operator $(x), belonging to the adjoint representation 
of SU(2). A gauge transformation g(x) which rotates 
$(a;) to (0,0, 1), or which diagonalizes <t(x) ■ a is called 
the abelian projection on &(x). g(x) can be singular in 
a configuration at the points where $(x) has zeros, and 
$(x) is not defined. 

The field strength F^, defined as pl| 

F^ = $(x) • 6^ - - (d„$(x) A DA(x)) • *(as) (2) 

is a colour singlet, and is invariant under non singular 
gauge transformations. 

In general eq. <Q) can be written as |2^ j 



r iu/ — O^Au d v A^ 

- - (c0(x) A d v $(x)) ■ $(x) 



(3) 



with 



A, = $<M£ 



(4) 



In the abelian projected gauge $(x) is constant, the sec- 
ond term in the right hand side of eq. (|3|) vanishes and 
the field F^ becomes an abelian field. 
Denoting by F* the usual dual tensor, 



and defining the magnetic current as 
i — 8" F* 

JfJ. — u r fly I 

it follows from eq.'s @, (§), (§) that 
d^J, = . 



(5) 



(6) 



(7) 



The magnetic charge is conserved, and defines a magnetic 
U(l) symmetry. 

The abelian projection g(x) can have singularities and 
as a consequence an additional field strength adds to the 
usual covariant gauge transform of G M „ pit] . 

After abelian projection 



G u 



gG^g' 



Gsniq 
flu 



(8) 



with G s ™ g 



$(x) (<9 M i;T 9 - dvA^A parallel to the 

colour direction $(x), and consisting of Dirac strings 
starting at the zeros of <&(x). The field configurations 
contain monopoles at the zeros of $(x), as sinks or 
sources of the regular field, and the strings carry away 
the corresponding magnetic flux. 

On a lattice (or in any other compact regularized de- 
scription in terms of parallel transport) the Dirac string 
reduces to an additional flux of 2n across a sequence of 
plaquettes, which is invisible [p3| . 

The mechanism relating confinement of colour to dual 
superconductivity of the vacuum advocates a sponta- 
neous breaking a la Higgs of the magnetic U(l) symmetry 
described by (|7j), which constrains the electric component 
of the field eq. (0) into flux tubes. 

All particles which have non zero electric charge with 
respect to the residual U{1) (0) will then be confined. 
There exist coloured states, e.g. the gluon oriented par- 
allel to $(x) which are not confined. This is a strong ar- 
gument that, if dual superconductivity is a mechanism of 
confinement at all, it must exist in many different abelian 
projections, as a manifestation of non abelian disorder. 

On the lattice the abelian gauge field corresponding to 
any given projection, or $(x), is extracted as follows [ p4[ . 

Let Ufj,(n) — g(n)Ufj,(n)g^(n) be the generic link after 
abelian projection. We adopt the usual notation U^n) = 

U^{n,t) = exp(iA^(n)), with A^n) = A^n) ■ a. 

The representation in terms of Euler angles has the 
form 



_ e «cv '3 e i7M' :r 2 e -''c v o-3 e i(/3 M + cv )<T3 



(9) 



e '^ 



•A* 



and 7 T is a vector perpendicular to the 3 axis. We assume 
the usual representation in which 173 is diagonal. 

For a plaquette, a similar decomposition can be per- 
formed, 



n 



(XV 



— e nw-<? e ^tT3 



(10) 

O^v = A M 0, y — /\ u 0^ up to terms 0(a 2 ). 6^ is the lattice 
analog of F^. The abelian magnetic flux is conserved 
by construction. A monopole appears whenever the flux 
entering five faces of a spatial cube adds to more than 2-k: 
then the flux through the sixth face is larger than 27r, but 
multiples of 2tt are invisible in the exponent. Formally 

M 



J l±v 



y flv 



27m 



y,v 1 



(ii) 



with — 7r < Q^y < 7r. A string through the sixth face 
takes care of the flux which has disappeared. 

We shall construct a disorder parameter for monopole 
condensation as the vev of an operator carrying non zero 
magnetic charge, /1. (/i) =/= will signal dual supercon- 
ductivity. 



III. THE DISORDER PARAMETER 

The disorder parameter will be constructed on the 
same lines as in ref.'s pjl2[|. 

An improvement exists with respect to ref . Mm , which 
consists in properly taking the compactness into account: 
in ref. |14| the approximation was that the field was 
treated as non compact. The same improvement was 
done in ref. [[12| with respect to ref. [11]. 

All the results presented in ref.'s [2C] already contain 
such improvement. 

We first analyze the case in which $(cc) is determined 
by the Wilson Polyakov line, i.e. the closed parallel trans- 
port to +00 along the time axis and back from —00 to 
the initial point via the periodic boundary conditions. 
For this choice, after abelian projection all the links 
Uo(n) along the temporal axis are diagonal, of the form 
U (n) =exp(iAl(n)a 3 ). 

Assuming for sake of definiteness the Wilson action we 
construct the operator /x(y, t) which creates a monopole 
at site y and time t with the following recipe (a similar 
construction can be made for other action). 

Let A M (x, y) be the vector potential describing the 
field value at site x of a static monopole sitting at y. We 
shall write it as 



?M 



A M (x,y)=Af(x,y)+VA(x,y), 



(12) 



with V- Af(x,y) =0. 

The first term describes the physical part of A M , the 
second term the classical gauge freedom. 

Let 1L0 be the electric field plaquette at time t. Then 
we define 



H = exp [-/3AS] , 



AS = \ E Tr Wo (".*)- 14 (",*)> 



Here 



U i0 (ft,t) = Ui(n,t)Uo(n + i,t) 

(i7 i (tJ,t+l)) t (J7o(rl,t)) t 



(13) 



(14) 



(15) 



is the electric field term of the action, and n' n is a mod 
ification of it, defined as 



(0 



n i0 (n,t)' = U l (n,t)U {n + it) (16) 

(Ul(fl,t+lrf(U (n,t)) 1 , 

U[{n,t + 1) = e jA ("^*("'*)- ff C/ 4 (n,t) (17) 

—iAYi(n,y)<b(n+i,t)-ff—iA(n+i,y)<h(n+i,t)-ff 

The disorder parameter is defined as (/i), or 



(m) 



/ iVU) e -^ s+As ^ 



IV>U)< 



-ps 



(18) 



It follows from the definition ([Tj) that adding AS to the 
action amounts to replace the term ITo at time t with 

The Ho (n, t) are the only terms in the action where 
the Uo{n,t) appear. In the path integral (jig) a change 

of variables U„(n,t) -> [/£(#,*) = U a (n,t)e iK( ^^ l - n ^- 3 
leaves the Haar measure invariant and reabsorbs the un- 
physical gauge factor of eq. ( |l7| ) , so that (fi) is indepen- 
dent, as it must be, of the choice of the classical gauge 
for the field produced by the monopole. 
Also a change of variables can be made 



Ui(ft, t + 1) -» Ui(n, t + i) e i^(fi,y)i(n+i,tyg 



(19) 



Again, this leaves the measure invariant, and brings 
Hio(n, t) to its original form. However in the plaquette 
Hij(ft, t+l) it produces the change Ui(n, t+l) — > f/j(n, i+ 
iyiAf t (ft,y)Mn+i,t)-s _ By the cons t m ction of Sect. |] this 
amounts to change, in the abelian projected gauge 



i(n,t + l) 



9ij{n,t + l) + 



(20) 



AiAfMffj-A^^y) 



or to add the magnetic field of a monopole. 

The same redefinition of variables reflects in the change 



IT (7M + l)^l4(n,i+l) 



(21) 



analogous to equation ([16]). Again the gauge factors 
e -iA(«,j/)#(»,t)-? ) e iA(ft+i,g)i(ft+i,tys are i rre i evant; s i n ce 

they can be reabsorbed in a redefinition of Uo(n,t + 1) . 
eIJ 4i f ,(™,j?)i>(n+M)-ff commutes with U (n+i, t+l), which is 
diagonal with it by definition of the Polyakov line abelian 
projection. 
In detail 

Il' l0 (n,t + 1) = U l {n,t + \)e lA ^ {n ^ {n+u) - s (22) 
U (rt + i,t + l) 
(C/ l (n,i + 2)) t (C/ (n,t + l)) t 
= Ui(n, t + l)U (n +i,t+ 1) 

iAYi(n,y)i(n+i,t)-S 

([/»(n,i + 2)) t (C/ (n,t + l)) t . 



A new change of variable can be done analogous to (19), 
exposing now a monopole at t+2 and producing a change 
Hio(n, t + 2) — » W i0 (n, t + 2). The procedure can be iter- 
ated. If an antimonopole is created at t + T, by an op- 
erator analogous to that of (|l3|), but with A^ — > — A± , 
then at time t + T the change cancels and the procedure 
stops. 

This shows that the correlation function 



indeed describes the creation of a monopole at y at time 
t and its propagation from t to t + T. This argument 
in this gauge is perfectly analogous to the argument for 
compact U(l) gauge theory [Q. The construction is the 
compact version of that of ref . Ej] . 
At large T, by cluster property 



V(T) ~ Aexp(-MT) + (/i) 2 



(24) 



where the equality (/_*) = (ft) has been used stemming 
from charge conjugation invariance. 

(/i) 7^ indicates spontaneous breaking of the U(l) 
magnetic symmetry defined in Sect. O eq. (pi), and hence 
dual superconductivity |20|. (/i) is the corresponding dis- 
order parameter. In the thermodynamic limit we expect 
(/i) ^ below the deconfining transition, (/i) = above 
it. At finite volume (fi) can not vanish for /3 > (3c with- 
out vanishing identically, since it is an entire function of 
(3. Only in the limit N s — » oo singularities develop p5| , 
and (fj,) can vanish. 

M is the lowest mass with quantum numbers of a 
monopole. In the Landau-Ginzburg model of supercon- 
ductivity it corresponds to the Higgs mass. When com- 
pared to the inverse penetration depth of the field, it can 
give information on the type of superconductor. We will 
discuss the determination of (fi) in the next section. 

For numerical reasons it will prove convenient to de- 
termine 



P=^log<M>, 



(25) 



which, by use of eq. (|1 
the two actions 



amounts to the difference of 



p=(S)s-(S + AS) s+ 



AS 



(26) 



V(T) = (fi(y,t + T)n(y,t)) 



(23) 



p contains all the informations we need. At finite temper- 
ature the lattice is asymmetric (Nf x Nt with Nt <C N s ), 
the quantities which can be computed are static and the 
vev of a single operator //, (/i) must be directly com- 
puted. Indeed there is no way of putting a monopole and 
an antimonopole at large distance along the t axis as we 
do at T = 0, since at T ~ T c , Nra is comparable to 
the correlation length. C*-periodic boundary conditions 
in time (U^n, Nt} = U*(n,0), where U* is the complex 
conjugate of U pi]) are needed. The magnetic charge 
is conserved. If we create a monopole say at t = 1, and 
we propagate it to t = Nt by the changes of variables 
described above, the magnetic charge at t = Nt will be 
different by one unit from that at t — 0, and this is in- 
consistent with periodic boundary conditions. With C*- 
periodic boundary conditions the magnetic field at Nt is 
opposite to the one at t = 0, since under complex con- 
jugation the term proportional to as in eq.(|9j) changes 
sign. By the change of variables eq. ( |19| ) a magnetic field 
in then added with opposite sign at Nt- This produces 
a dislocation with magnetic charge —1 at the boundary 
which plays the role of the antimonopole in eq. (E3h . 



With a generic choice of the abelian projection different 
from the Polyakov line, we can define the operator p in 
a similar way, by sending 



n i0 (n,i)^l4(n,i) , 



(27) 



according to eq. (|16[). 

Again to demonstrate that a monopole is created at 
t + 1 we can perform the change of variables, eq. (19), 



and expose a change of the abelian magnetic field at t + 1 
given by eq. (f20|). 

However now the resulting change of Iloi (n, t + 1) is 

I4(»f,t+1) = ^(fM+l)e M " (fy7) * ( " +M) '* (28) 

C/ (n + z,i+l) 

([/,(«,£ + 2)) t ([/ (n,i + l)) t 
= £/i(n, i + l)L/ (n + U + 1) 

((C/ (n + t, t + l))y e ^ h A(n,y)*(n+tt)-a 

U (n + i,t+l) 

(t/. i (n,i + 2)) t (C/ (n,t + l)) t . 
The change of Ui(n, t + 2) is by a factor on the right 

(!7o(n + M+ i))t e *A£(fl,fl*(fl+M)-* (29) 

[7 (n + i,t + f) . 

C/o does not commute with e iA ^ i ( n '^*( n+l '*)' cr , as it was 
in the case of $(n, t) in the direction of the Polyakov line. 
This looks at first sight as a complication, but it is not. 
Indeed the abelian projected phase of a product of links 
is the sum of the abelian projected phases of the factors, 
to 0{a 2 ). From eq. @ it follows 0{a 2 ) 



-~T1 -* - 1 -~T2 -. - 2 



jf T -a J(ll+"tl)<?3 



(30) 



Hence the abelian phases of Uq, Uq in ( |29| ) cancel 0(a 2 ) 
and the abelian projected field of the modified plaquettc 
at time t + 2 is again changed according to eq. (p0|). 



IV. NUMERICAL RESULTS FOR p 

We will determine the temperature dependence of p on 
an asymmetric lattice N% x N t (N s » Nt). 

For reasons which will be clear in what follows we 
will distinguish between abelian projections in which the 
operator $(n, t) which defines the monopoles is explic- 
itly known, and projections (like the so called Maximal 
Abelian) in which the projection is fixed by a maximizing 
procedure, and <&{n,t) is not explicitly known. 

In the first category we studied the following projec- 
tions. We will define the operator <$>(x) = $(x)/|$(x)| 
starting from an operator O which is an element of the 
group, by the formula 

O = Co + i${x) ■ a 



• O is connected to the Polyakov line L(n,t) = 
lfc 1 fo(ri,* / )l4'= 1 o c/ b(rf,* / ) as follows^ 

0(fl,t)=U^ 1 Uo(n,t')L*(ft,0)II t t z} U (n,t') ; (31) 

• O is an open plaquette, i.e. a parallel transport on 
an elementary square of the lattice 

G(n) = Uij(n) = U.inpjin + i) (32) 

(Uiin + mHUjin)? ; 

• "butterfly" projection, where the projecting oper- 
ator is 

0{n) = F(n) = U x (n)U y (n + x) (U x (n + y)) f (33) 
(Uy(n))^ U z (n)U t (n + z) 

(U z (n + i)) f (Ut(n)^ . 

The trace of F is the density of topological charge. The 
projection defined in (p^) is the Polyakov projection on 
a C*-periodic lattice. 

From eq. (p5J) and the condition (/i(/3 = 0)} = 1 we 
obtain 



</z(/3)) = exp / p((3')d(3' 



(34) 



If (//) (defined in any abelian projection) is a disorder 
parameter for the deconfining phase transition, we expect 
that in the thermodynamic limit (N s — > oo, Nt constant) 
p goes to a finite bounded value in the strong coupling 
region, i.e. in the region below the deconfining transition. 
In the weak coupling region (p) should go to zero in the 
same limit, i.e. p must go to — oo. In the critical region 
we expect an abrupt decrease of (//) , and hence a negative 
sharp peak in p. 

A few details about numerical computation. Accord- 
ing to eq. (Eq), p is the difference between two actions: 
the standard SU(2) Wilson action and the "monopole" 
action S + AS. 

For the Wilson term simulation can be performed by 
using an heat-bath algorithm. This is not possible in 
the case of the "monopole" action. Consider for example 
the Polyakov projection and a single monopole operator 
p(y,0). In the updating procedure, we can distinguish 
the following four cases: 

1. update of a spatial link at t ^ 0, 1. The plaquettes 
involved have Wilson's form and the variation of 
the "monopole" action is linear with respect to the 
link we are updating; 



by * we indicate the complex conjugation operation. 



2. update of a spatial link at t = 0, 1. Although some 
plaquettes are modified by the monopole term, the 
variation of the modified action S + AS is again 
linear with respect to the link, because the field <& 
does not depend on the link we are changing; 

3. update of a temporal link at t ^ 0. The local vari- 
ation of the action is linear, but the change also 
induces a change of the Polyakov loop, i.e. of $, 
according to eq. (pl|), so that there is an effect on 
the action which is non linear; 

4. update of a temporal link at t = 0. We can not de- 
fine a force, because due to the change of the corre- 
sponding Polyakov loop the change of the modified 
action is non-linear. 

In order to perform numerical simulations in the case of 
the system described by the "monopole" action, an ap- 
propriate algorithm could be Metropolis; however, this 
method can have long correlation times. In view of im- 
proving decorrelation, we have performed simulations by 
using an heat-bath algorithm for the update of the spa- 
tial links and a Metropolis algorithm for the update of 
the temporal links. 

Similar techniques can be used for the other projec- 
tions we have investigated: in all cases, we have chosen 
to use the heat-bath updating when the contribution to 
the action is linear with respect to the link we are chang- 
ing and the Metropolis algorithm when it is not. As a 
test we verified that the mixed update correctly works 
for the Wilson action. 

The simulation was done on a 128-node APE Quadrics 
Machine. We used an overrelaxed heat-bath algorithm to 
compute the Wilson term of eq. d29) , and a mixed algo- 
rithm as described above for the other term. Far from 
the critical region at each (3 typically 4000 termalized 
configurations were produced, each of them taken after 4 
sweeps. The errors are computed with a Jackknife anal- 
ysis to the data binned in bunches of different length. As 
error we took the maximum of the standard deviation as 
a function of the bin length at plateau. In the critical re- 
gion a higher statistic is required. Typically the Wilson 
term is more noisy. Thermalization was checked by moni- 
toring the action density and the probability distribution 
of the trace of the Polyakov loop. 

The discussion of Sect. Ill implies that different choices 
for A A f are equivalent: eq. ( po| ) shows that only the 
magnetic field of the monopole determines the value of 
(p) . In our simulation we used the Wu-Yang form of A± . 
We checked that the Dirac form (with different position 
of the string) gives compatible results. 

In simulations of S + AS we found that correlation 
times are small and under control for Nt = 4. For N t — 6 
in the critical region thermalization problems arise and 
modes with long correlation time appear. For this reason, 
we have used mainly lattices with N t — 4. 

Fig. PP shows the typical behaviour of p for different 
abelian projections, for a lattice 12 3 x 4. The negative 



peak occurs at the expected transition point, (3c P7[ . 
Below 0c the different projections are indistinguishable 
within errors, suggesting that different monopoles behave 
in the same way. 

Fig. shows the comparison with a 18 3 x 6 lattice. 
The peak is displaced at the correct (3c, showing that it 
is not an artifact but it is related to deconfinement of 
colour @. 
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Since different projections give indistinguishable re- 
sults, for sake of simplicity we shall only display the pla- 
quette projection in the following figures. 

Fig. g shows the dependence of p on N s at fixed Nt — 
4. The qualitative behaviour does not change when we 
increase the lattice size. We now try to understand the 
thermodynamic limit |q,|l2|. 

In the strong coupling region (cfr. fig. 0) p seems to 
converge to a finite value, which is consistent with at 
low /?'s. Eq. (p4T) then implies that (/z) ^ in the infinite 
volume limit in the confined phase. 

The weak coupling region is perturbative. An estimate 
of p is the minimum on the ensemble of the configurations 



U and is given by the action of classical solutions of the 
system described by S + AS: 



p — > minis'} — minis' + AS} 

/3— >oo L U U 

= - min|S + AS} , 

u 



(35) 



since minjy|S} = 0. 

In other systems, where the same shifting procedure 
has been applied and studied, this asymptotic value has 
been analytically calculated in perturbation theory with 
the result [§JI§ 



-cN s 



(36) 



where c and d are constants, i.e p goes linearly with the 
spatial dimension. 

In SU(2) we are unable to perform the same calcula- 
tion and we have evaluated the nrinimum min;y{S-|- AS} 
numerically. Some technical remarks on the numerical 
procedure. An oversimplified strategy would be to start 
from a random configuration and then decrease the ac- 
tion by Metropolis-like steps in which the new configura- 
tion is accepted only if its action is lower. However this 
procedure will not work, because of the presence of local 
minima where often the procedure stops. A way to over- 
come this difficulty is to perform an usual Monte-Carlo 
simulation where (3 is increased indefinitively during the 
simulation p8,M. This is equivalent to freeze the sys- 
tem. 

We found useful to integrate the two strategies. Firstly 
we freeze the system increasing (3 in the following way: 

1. we thermalize the system at a reasonable (e.g. 

= 10); 

2. we increase (3 by a fraction 1/200 and at the new 
value of j3 we perform a number of sweeps (typically 
200), looking for the corresponding minimum of the 
action; 

3. we iterate the step 2 until the minimum of the ac- 
tion looks stable along a larger number of sweeps 
(typically 5000). 

When this procedure becomes inefficient (typically for 
(3 rs 10 6 ), we go to a Metropolis- like minimization, which 
is stopped when the action stays constant within errors. 

The result is shown in fig. || for the plaquette projec- 
tion. It is consistent with the linear dependence of eq. 
( Pq ) with c ~ 0.6 and d ~ —12. Thus in the weak cou- 
pling region in the thermodynanric limit p goes to — oo 
linearly with the spatial lattice size and 



(/'•) 



N, 



Ae (-cN s +d)0 _> Q) p >l3i 



■c 



(37) 
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FIG. 5. p vs. N s (N t = 4) at (3 — oo in the plaquette 
projection. Data are obtained by numerical minimization of 

S + AS. 
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To sum up, (/i) is different from zero at least in a wide 
range of (3 below (3c and goes to zero exponentially with 
the lattice size for (3 > (3c- The strong coupling region 
and the weak coupling one must be connected by a de- 
crease of (n) ; the sharp peak of p signals that this decline 
is abrupt and takes place in the critical region. 

To understand the behaviour of p near the critical 
point we shall use finite size analysis. By dimensional 
argument 



(u) = N~ s / V $ f -1. - ^ 



(38) 



where a and £ are respectively the lattice spacing and the 
correlation length of the system. 
Near the critical point, for /3 < (3c 



t<x{pc-0)~ 



(39) 



where v is the corresponding critical exponent. In the 
limit N s ^> Nt and for a/£ <C 1, i.e. sufficiently close to 
the critical point we obtain 



( Ji ) = Nr s/v *(Ny*'(pc-0),o,o 



or equivalently 



N, 



l/v 



f (Nl' v {(3 c - (3) 



(40) 



(41) 



The ratio p/N s is a universal function of the scaling 
variable 



x = N l J v ((3c - (3) 



(42) 



Since critical values of (3 and critical indices of SU(2) 
pure gauge theory are well known |]27[ , we can check how 
well scaling is obeyed by plotting p/N s ' 1 ' as a function 
of x. 

Fig. ra shows the quality of the scaling in the plaque- 
tte projection for (3c — 2.2986 and v = 0.63. Similar 
qualitative results have been obtained for the Polyakov 
projection. 

As a further check, we can vary v and try to estimate 
"by eye" sensitivity of our data to this exponent. We 
obtain that in both projections the scaling relation is 
satisfied within errors for 0.57 <v< 0.67. 

In the thermodynamic limit in some region of (3 below 
the critical point we expect 



(H) ex ((3 C - (3) d 



which implies 



N 



l/v 



X 



(43) 



(44) 



Using eq. (44) it should be possible in principle to de- 
termine v, 5 and (3c- Our statistic is not enough accurate 
to perform such a fit. However, we can determine 8 using 
as an input (3c , v, which are known, by parameterizing 
p in a wide range by the form 



P 



N, 



l/v 



X 



(45) 



where c is a constant, as suggested by fig. |6J. 

Our best fit| gives 8 = 0.24 ± 0.07 in the plaquette 
gauge and 8 = 0.12±0.04 in the Polyakov gauge. The 
reduced x 2 is order 1. 



"Fits have been performed by using the Minuit routines. 



This concludes our argument about the thermody- 
namic limit (N s — > oo). The deconfining phase transition 
can be seen from a dual point of view as the transition 
of the vacuum from the dual superconductivity phase to 
the dual ordinary phase. That feature seems to be inde- 
pendent of the abclian projection chosen. 



V. THE MAXIMAL ABELIAN PROJECTION 

There are abelian projections which are not explicitly 
defined by an operator <I>, but by some extremization 
procedure. The prototype is the Maximal Abelian pro- 
jection, which is defined by maximizing numerically the 
quantity 



S l u(W)=^tr[[/ /J (n)<T 3 t/ / i(n)<73] 



(46) 



with respect to gauge transformations [|18|,[50| . 

The Maximal Abelian projection is very popular since, 
in the projected gauge all links arc oriented in the abelian 
direction within 15%, and therefore all observables are 
dominated by the abelian part within 85%. This fact 
is known as abelian dominance |3l| , and could indi- 
cate that the abelian degrees of freedom in this projec- 
tion are the relevant dynamical variables at large dis- 
tances. Moreover, out of the abelian projected configura- 
tions, monopoles seem to dominate observable quantities 
(monopole dominance [p2[). 

With our approach we have a technical difficulty to 
determine p via S + AS (eq. (|26|)). At each updating 
the operator <I> and S + AS are only known after max- 
imization. Accepting or rejecting an updating therefore 
requires a maximization, and the procedure takes an ex- 
tremely long computer time. 
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Therefore in order to study this abelian projection we 
have to explore the possibility of measuring (fi) directly, 
and to confront with the huge fluctuations coming from 
the fact that \x is the exponential of a sum on a space 

N 3 / 2 



volume and typically fluctuates as 
the following strategy [Of 



We adopt 



1. we study the probability distribution of the quan- 
tity 



log/i = -/3(AS) ; 



(47) 



2. we reconstruct (fi) from the log/i distribution by 
means of cumulant expansion formula truncated at 
some order. 

This procedure should be compared with that of ref. |33| . 
If we have a stochastic variable X distributed with 
probability p(X) we have 



dXe m- W)p(X) = e^>= 



-c„ 



(48) 



(X) is the mean value and C n is the n-th cluster. For 
example, if we call A = X — (X) 



d =0 : 

C2 = (A 2 ) ; 

C 3 = (A 3 ) ; 

C7 4 = (A 4 )-3(A 2 )(A 2 ); 

C 5 = (A 5 ) - 10(A 3 )(A 2 ) 



(49) 



If log/i were gaussianly distributed with mean value m 
and standard deviation a, we would have 



(/i) = exp ( to + 



(50) 



In order to check the method we have explored the 
cluster expansion for the logp distribution in the projec- 
tion we have already studied by means of the quantity p. 
Fig. |7] shows a comparison between \og(p), taken from the 
integration of p data, and cluster expansions truncated 
at different orders. The first and the second cluster are 
insufficient to account for the right behavior of log(/i), 
whereas with the third cluster added the two determina- 
tions are consistent. Moreover the fourth cluster is zero 
within statistical errors. It seems that one can estimate 
(p) with a cluster expansion truncated at the third order. 
As a rule, the higher clusters are quite noisy and error 
bars grow with increasing order. Therefore this kind of 
estimation requires a very high statistics. For this rea- 
son numerical determination of log{p) in the Maximal 
Abelian projection is possible, but very time consuming. 
Our data are displayed in fig. g| showing that monopoles 
in the Maximal Abelian projection behave in the same 
way as monopoles in other projections. 

For this kind of simulations, we have used a standard 
overrelaxed heat-bath algorithm. For each value of (3 
we performed about 50000 measurements, each of them 
taken after 8 sweeps. In order to improve the statistics, 
we have considered eight symmetric different position of 
the monopole (namely we have inserted the monopole at 
the center of each optant of a cartesian coordinate sys- 
tem with the origin at the center of the lattice); data 
corresponding to each position are analyzed separately 
with the method exposed in the previous section and our 
best value is the weighted average of the eight measure- 
ments. Putting more monopoles would not improve the 
statistics, since strong correlations appear whenever the 
distance is shorter than the correlation length. Also these 
simulations have been performed on a 128-node APE 
QUADPJCS machine. 



VI. CONCLUSIONS 



when T — > Tc can be determined. The first two quan- 
tities are known independently and our determination is 
consistent with others. As for 6, defined by 



^ tZt c V T c 



(51) 



it is 0.20 ±0.08. Different abelian projections (plaquette, 
Polyakov, "butterfly") give results which agree with each 
other. 

Our technique proves difficult for the Maximal Abelian 
projection, but a direct determination of {p) looks con- 
sistent with other projections. 

In conclusions 

1. Dual superconductivity is at work in the confined 
phase, and disappears at the deconfinement phase 
transition. 

2. This statement is independent of the abelian pro- 
jection defining the monopoles. 

Further theoretical effort is needed to understand the 
real symmetry breaking in the deconfined phase, or in 
the dual description of QCD. 

Similar results for SU(3) will be presented in the com- 
panion paper. 

Finally we stress that, whatever topological excitations 
are responsible for colour confinement, counting them is 
not a right criterion to detect disorder. Only the vev of 
an operator carrying the appropriate topological charge 
can be a legitimate disorder parameter. 
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We have constructed a disorder parameter (p) detect- 
ing condensation of monopoles of non abelian gauge theo- 
ries defined by different abelian projections. The param- 
eter is the vev of an operator which creates a magnetic 
charge, (p) =/= signals dual superconductivity. The 
same construction has been tested in many known sys- 
terns @g|g|. 

We measure by numerical simulations (p), or better 

p = -]— log(/i), which contains all the relevant informa- 
tions and less severe fluctuations. 

An extrapolation to thermodynamic limit (infinite spa- 
tial volume) is possible. 

The system behaves as a dual superconductor in the 
confined phase, and has a transition to normal at the 
deconfining phase transition, where (p) — » 0. 

The deconfining (3c and the critical index v as well as 
the critical index S describing the way in which (p) — ► 
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